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BACKGROUND OF THE INVENTION 



1. Related Application 

This application claims priority from co-pending provisional application Serial No. 
60/152,391 filed September 3, 1999. 

5 2. Field of the Invention 

This invention relates generally to the field of fluid dynamics. More particularly, the 
invention is a method for accurately and efficiently analyzing anisotropic turbulent flows of both 
gases and liquids. 

3. BACKGROUND 

The study and engineering application of hydraulics or fluid mechanics dates to the dav^n of 
civilization. The first irrigation canals were constructed before 5,000 BC. Theoretical analysis of 
fluid motions began in the modem scientific sense with Newton's Principia (1687) and his laws of 
motion. Bernoulli (1738) established a relationship between the pressure and velocity of an ideal 
fluid which became his famous theorem. d'Alembert (1744) deduced that the steady flow of an 
ideal fluid about a body produces no drag (d'Alembert' s Paradox). Euler (1755) founded modem 
fluid mechanic analyses by deriving the differential equations of motion (continuity and momentum 
conservation) for an ideal (inviscid) isotropic fluid. Navier (1822) and Stokes (1845) 
independently derived the viscous terms necessary to extend Euler' s equations to isotropic, viscous 
fluids. The equation of state and the conservation of energy equation (First Law of 
Thermodynamics) written utilizing Fourier's (1822) law of thermal energy transfer completed the 
classical Navier Stokes equations governing isotropic fluid motions. 

Reynolds (1883) experimentally showed that there are two modes of fluid motion, laminar 
and turbulent motion, which occurs at large values of the dimensionless Reynolds number. 
Reynolds (1894), by considering a time average motion and a time dependent turbulent motion, 
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showed that turbulent motions introduce additional turbulent (Reynolds) stresses into the Navier 
Stokes equations, greatly increasing the number of flow properties required to describe turbulent 
flows. 

Despite over a century of intense aircraft and aerospace developments, wind tunnel and 
5 laboratory experiments, theoretical analyses and numerical computations based on the Navier 

Stokes equations, the determination of a general theory of turbulent motion remains one of the last 
great unsolved problems of classical physics. 

Maxwell (1858) founded modem kinetic theory when he introduced his new concept of the 
molecule velocity distribution function and deduced its equihbrium form. Maxwell (1867) derived 
□ 1 0 the Maxwell transport equations and showed that the coUision change integrals could be solved 
n) analytically without knowing the molecule velocity distribution function for Maxwell molecules. 

£ Boltzmann (1872) derived his integro-differential equation for Maxwell's molecular velocity 

U distribution function and solved for his famous H-theorem, which proved that Maxwell's 

J. equilibrium form was correct. Chapman (1916), using Maxwell's transport equations, determined 

1 5 accurate general formulae for the gas transport coefficients. Enskog (1917) gave a general solution 
method for he Boltzmann equation, which showed that the Euler equations were the first and the 
Navier Stokes equations were the second approximate solutions to the Boltzmann equation and 
gave gas transport coefficients identical to Chapman's. Both the Chapman and Enskog analyses 
were for isotropic perfect gases, being based on perturbations of Maxwell's isotropic equilibrium 
20 solution. Their independent analyses cemented the belief in both the scientific and particularly the 
engineering professions that the Navier Stokes equations were generally and universally valid and 
that the correct way of analyzing both laminar and turbulent flows was through the Navier Stokes 
equations. This belief has been engraved in stone in both the scientific and engineering literature 
throughout the 20th century. 
25 Burnett (1935) derived the third approximate solution to the Boltzmann equations using 

Enskog' s solution method. Grad (1949) derived a thirteen moment method of solving 
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Boltzmann's equation which closely followed Enskog's methods. Both of these analyses have 
been proven to lack generality and have not been significant advances of the Navier Stokes theory. 
Bird (1963) introduced the Direct Simulation Monte Carlo (DSMC) numerical method for 
calculating rarified gas flows. Yen (1966) showed that the directional thermal energies 
(temperatures) were vastly different in shock waves and that the longitudinal temperature overshot 



its downstream value for shock Mach numbers greater than J— (1.34). Elliot and Baganoff 



(1974) showed that the Navier Stokes normal stress relationship and the Fourier energy flux 
component ratio were valid only in sound waves and were invalid at the shock end points. Elliot 

(1975) showed that the Navier Stokes normal stress relationship was invalid everywhere in shock 
waves and incompatible with the directional thermal kinetic energy moments of the Boltzmann 
equation. DSMC numerical shock calculations confirmed all of these analytical predictions and 
numerically illustrated the anisotropic not isotropic fluid features of many simple gas flows. 

Kliegel (1990), following a suggestion of Maxwell (1867) that near equilibrium flows were 
anisotropic and should be represented by an anisotropic, not isotropic Maxwellian, performed the 
first general anisotropic solution analysis of the Boltzmann equation. Kliegel (1990) showed that 
the Euler-Navier Stokes-Bumett equation sequence was not the correct approximate solution 
sequence for the Boltzmann equation. He showed that the correct gas dynamic equation set was an 
anisotropic fluid seven equation set for the density, three fluid velocity components and three 
directional thermal kinetic energies (temperatures), not the classic isotropic fluid five equation 
Euler-Navier Stokes set for the density, three fluid velocity components and total temperature. He 
also gave shear and directional energy flux relationship having the correct Mach number 
dependence for both sound waves and weak shocks. Kliegel' s anisotropic equations resolved 
d'Alembert's paradox giving a profile pressure drag associated with pushing an ideal (shearless, 
energy fluxless) fluid about a body. They also correctly predicted the directional thermal energy 
separations and overshoots and the thermal energy flux component changes occurring in shock 
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waves as predicted by Yen (1966) and numerical DSMC shock calculations. Bird (1994) 
summarizes recent DSMC calculational capabilities. Chen and Jaw (1998) presented a recent 
summary of classical isotropic fluid Navier Stokes based turbulent flow modelling. 

The present disclosure correctly consolidates all previous fluid dynamic and kinetic theory 
5 analyses and extends it to define a new method of analyzing and computing anisotropic turbulent 
flows. 

The prior art cited above may be found in the following references: 

1. 1687, Newton, I. Philosophiae Naturalis Principia Mathematical Oxford. 

2. 1738, Bernoulli, D., Hydrodynamica sive de Viribus et Motibus Fluidorum Commentarii, 
1 0 Beriin. 

3. 1744, d' Alembert, J., Traite de Vequilibre et du mouvement desfluides pour servir de suite au 

traite de dynamique, Paris. 

4. 1755, Euler, L., Principles genereaux du mouvement desfluides, Histoire de 1' Academic de 

Berlin. 

15 5. 1822, Navier, L.M.H., Memoires de V Academic des Sciences de VInstitut de France, V 6, p. 

389, Paris. 

6. 1822, Fourier, J., Theorie analytique de la chaleur, Paris. 

7. 1845, Stokes, G.G., On the Theory of the Internal Friction of Fluids in Motion, Transactions of 

the Cambridge Philosophical Society, V8, p. 287, Cambridge. 
20 8. 1858, Maxwell, J.C., Illustrations of the Dynamical Theory of Gases, Philosophical Magazine, 

V 19, p.l9andV20,p.21. 

9. 1867, Maxwell, J.C., On the Dynamical Theory of Gases, Philosophical Transactions of the 

Royal Society, V 157, p. 49. 

10. 1872, Boltzmann, L., Weitere Studien Uber das Wdrmegleichgewicht unter Gasmolekiilen, 
25 Styungsberichte, Akad. Wiss, Vienna, Part U, V 66, p. 275. 
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1 1. 1883, Reynolds, O., An Experimental Investigation of the Circumstances which Determine 

whether the Motion of Water Shall Be Direct or Sinuous, and the Law of Resistance in 
Parallel Channels, Philosophical Transactions of the Royal Society, V174, p. 935. 

12. 1895, Reynolds, O., On the Dynamical Theory of Incompressible Viscous Fluids and the 

Determination of the Criterion, Philosophical Transactions of the Royal Society, V186A, p. 
123. 

13. 1904, Prandth, L., Uber FlUssigkeitsbewegung bei sehr kleiner Reibung, Proceedings of the 

Third International Mathematical Conference, Heidelberg, Leipzig. 

14. 1916, Chapman, S., The Kinetic Theory of Simple and Composite Monotonic Gases: Viscosity, 

Thermal Conduction, and Diffusion, Proceedings of the Royal Society London, A, V 93, 

p.i. 

15. 1917, Enskog, D., Kinetische Theory der Vorgdnge in mdssig verdunnten Gasen, Almquist & 

Wilsells, Uppsala, Sweden. 

16. 1935 Burnett, D., The Distribution of Molecular Velocities and the Mean Motion in a Non- 

Uniform Gas, Proceedings London Mathematical Society, V 40,p. 382. 

17. 1949, Grad, H., On the Kinetic Theory of Rarefied Gases, Convention of Pure and Applied 

Mathematics, V 2, p. 33 1 . 

18. 1963, Bird, G.A., Approach to Translational Equilibrium in a Rigid Sphere Gas, Physics of 

Fluids, V 6, p. 1518. 

19. 1966, Yen, S-M, Temperature Overshoot in Shock Waves, Physics of Fluids, V9, p. 1417. 

20. 1970, Chapman, S., Cowling, T.G., Mathematical Theory of Non-Uniform Gases, Third 

Edition, Cambridge University Press. 

21. 1974, Elliot, J.P., Baganoff, D., Solution of the Boltzmann Equation at the Upstream and 

Downstream Singular Points in a Shock Wave, Journal of Huid Mechanics, V 65, p. 603. 

22. 1975, Elliot, J.P,, On the Validity of the Navier-Stokes Relation in a Shock Wave, Canadian 

Journal of Physics, V 53, p. 583. 
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23. 1990, Kliegel, J.R., Maxwell Boltzmann Gas Dynamics, Proceedings of the 17* International 

Symposium on Rarefied Gas Dynamics, Aachen 1990, Edited by Alfred E. Beylich, VCH, 
New York. 

24. 1994, Bird, G.A., Molecular Gas Dynamics and the Direct Simulation of Gas Flows, Oxford 

University Press, New York. 

25. 1998, Chen, C-J. Jaw, S-Y, Fundamentals of Turbulence Modeling, Taylor & Francis, 

Washington. 
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SUMMARY OF THE INVENTION 

This disclosure concerns a simple method of deriving the correct kinetic anisotropic fluid 
turbulent flow equations from the Boltzmann equation and the correct thermal and turbulent 
moment closures to yield a new, closed set of fluid dynamic equations describing anisotropic 
5 turbulent motions, solving the century old quest for the solution to this classical problem. 

Current turbulent flow analyses are based on the isotropic Navier Stokes equations, turbulent 
fluctuation moments of these equations and modeling of the unknown turbulent moments, identifying 
introduced unknown modeling coefficients by experiment. This results in a variety of semi-empirical 
turbulent flow analyses based on different Navier Stokes moment relationships satisfying mass, 
1 0 momentum and total energy conservation for an isotropic fluid, calibrated to approximately match a 
limited number of experiments. The limitations of current Navier Stokes based turbulent analyses are 
well known and appear to be irresolvable after a century of intense analytical and experimental effort. 

The present invention provides a general, closed, anisotropic kinetic turbulence theory for gases 
and liquids based on new solutions of the Maxwell moment equations of the Boltzmann equations. 

1 5 These solutions provide a closed initial equation set for the four time average fluid mechanic variables, 
the sixteen time average thermal motion correlation and the sixteen time average turbulent motion 
correlations listed in Table I. Higher order closure sets involving higher order (fourth and fifth, etc.) 
thermal and turbulent moment sets are easily generated allowing solutions of turbulent flow problems 
to higher and higher accuracy if desired or needed for special problems. These higher order sets will 

20 not be further discussed since they add only unnecessary detail to the disclosure and reveal Uttle that is 
new. 

An algebraic differential equation solution method described in my co-pending application 
entitled "Method for Algebraically Solving Differential Equations, including Stiff Equations, to High 
Accuracy" may be used to reduce the thirty six governing differential equations to the solution of as 
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few as six differential equations (continuity, momentum, total energy and turbulent energy) plus thirty 
algebraic relationships. The actual equation set solved depends on the problem and accuracy desired. 

The presented analytic anisotropic turbulence theory is mathematically correct, involving no 
semi-empirical modeling. The even order higher order moment closures involved are general and 
5 correct to known accuracy. The analysis can be used to solve the infinite set of Maxwell moment 
equations to any level desired to verify the accuracy of the initial closure set. 

The anisotropic kinetic turbulence theory presented allows the calculation of turbulent motions 
with the same accuracy and computer resources as current laminar flow calculations. 

The new teachings of the present invention may be summarized as follows: 

10 1. A closed turbulence theory may be simply derived directly from the Boltzmann 



equation without turbulent correlation modelling or experimentally determined 



modelling coefficients. 



2. 



Turbulent flows are totally in viscid flows, unchanged by molecular collisions. 



3. 



Turbulence is universal, having the same governing equations and the same physical 



15 



structure in both gases and liquids because turbulence is an inviscid flow 



phenomena described entirely by inviscid convective flow equations, independent of 



collisions and molecular interaction laws. 



4. 



The fact that time average turbulent flow correlation between turbulent flow 



fluctuations and density or thermal velocity correlations are zero because there is no 



20 



physical (coUisional) connection between these fluctuations, since both the turbulent 



flow Auctions and density are conserved in collisions. 



5. 



A method of analyzing the Maxwell moment equations to obtain a closed set of 



equations which describe time average turbulent flows to any desired accuracy level. 
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• 



All higher moment equation above the closure level can be solved in even/odd sets of 
two, thus yielding solutions to all Maxwell moment equations with error of Burnett 
order or less. 

6. Direct predictions of turbulent transition locations, turbulent energy growth rates, 
shears and other thermal and turbulent flow quantities of physical importance can 
now be derived or calculated from solution of the presented equations with the same 
accuracy and computer resources as current laminar flow analyses or calculations. 

7. Universal even order thermal and turbulent moment closure relationships exist and 
are explicitly given, allowing the termination and solution of the Maxwell moment 
equations at any desired odd moment equation set, and solution of the Boltzmann 
equation to any desired approximation. 
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DETAILED DESCRIPTION OF THE INVENTION 

In the following description, for purposes of explanation and not limitation, specific details 
are set forth in order to provide a thorough understanding of the present invention. Hov^ever, it will 
be apparent to one skilled in the art that the present invention may be practiced in other 
5 embodiments that depart from these specific details. In other instances, detailed descriptions of 
well-known methods and devices are omitted so as to not obscure the description of the present 
invention with unnecessary detail. 

The Boltzmann equation may be written in fluid dynamic collision rate form: [See 
Appendix, Equation set 1]; 

1 0 where = [pF(v)] is the Lagrangian convective motion derivative, p is the fluid density, o) is the 
molecular collision rate, F(v) is the molecular velocity (v) distribution function and F(v)* is the 
molecular velocity distribution function collided into F(v). There is no approximation in writing the 
Boltzmann collision integral in this format, since the required thermal velocity moments of F(v)* 
are correctly calculated from the full Boltzmann coUision integral (Kliegel (1990)). The molecular 

1 5 collision rate is 

for Maxwell molecules, where K is a molecular collision constant. 

The illustrative derivation will be done for Maxwell molecules for simplicity, since the 
colUsion change integrals are known exactly for Maxwell molecules. The same derivation can be 
performed for any molecular model, it is just more algebraically complex and adds nothing to the 
20 teachings. 

Both the original Enskog (1917) solution as a perturbation about the equilibrium 
Maxwellian and the Kliegel (1990) solution as a perturbation about a non-equilibrium anisotropic 
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Maxwellian are Navier Stokes order accurate Burnett order error terms) and do not depend on 

density gradients. Thus, density gradient corrections to the molecular velocity distribution function 
and its moments belong to the Burnett order of corrections. This important fact, although present in 
the result of both analyses, is expUcitly stated by neither author nor utilized to further extend or 
5 simplify the analyses. 

Multiplying the Boltzmann equation by the molecular velocity product v|v2 V3, where /, m, 
n are integers, and integrating over all possible velocities (indicated by the bar) yields the Maxwell 
molecular velocity moment equation: [See Appendix, Equation Set 2]. Derivation of the 
anisotropic turbulent flow equations from this equation will be discussed later. 

1 0 Consider the molecular velocity (v) to be composed of a mass mean fluid flow velocity (u) 

and a thermal motion velocity (c), then the Cartesian components of these velocities are related by 

Vi=Ui-hCj 

Multiplying the Boltzmann equation by the thermal velocity product cjc^c", where /, m and n are 
integers, and integrating over all possible thermal velocities (indicated by the bar) yields the 
1 5 Maxwell thermal velocity moment equation: [See Appendix, Equation Set 3], 

For odd moment (/ + /w + « = 2N -h 1, N integer) equations, the coefficients of the density 
gradient terms are even moments and products of lower even moments, which have equilibrium 
values. Thus these even moment coefficients must be zero to Burnett order for the density gradient 
terms to be of Bumett order. This requirement yields the general even order (fourth, sixth and 
20 higher) density gradient independent thermal moment closure relationships: [See Appendix, 
Equation Set 4]. 

The fourth order thermal moment closure relationships can also obtained from Kliegel 
(1990) by simple calculation but he did not state these results nor utilize them in any manner. 
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Substituting these relationships into Equation Set 3, one obtains the closed odd order (third, fifth 
and higher) density gradient independent thermal moment equation. [See Appendix, Equation Set 



5 Equation Set 6]. 

After using the fourth order moment density gradient independent closure relationships: 
[See Appendix, Equation Set 7], the closed directional thermal energy flux and structure equations 
become: [See Appendix, Equation Set 8] where the fourth order thermal moments have been 
replaced by products of second order moments through the fourth order density gradient 
1 0 independent closure relationships. 

The above set of equations are the correct anisotropic gas dynamic flow equations with 
neglected terms of Burnett order in the convective terms of the closed thermal energy flux and 
structure equations. All other equations in the set are exact. These equations resolve the known 
deficiencies of the Euler-Navier Stokes-Bumett equation sets, d' Alembert's paradox, the normal 
1 5 stress incompatibility, the incorrect Mach number dependency of the Navier Stokes normal stress 
and Fourier energy flux components, and the incorrect coordinate rotation variance of the Burnett 
equations. These equations are valid for both sound waves and weak shocks, and thus will 
accurately predict all flows except in those regions of high chcinge rates in distances comparable to 
the mean free path length, such as occur in strong shock waves and at boundaries in very rarified 



The twenty moment equations may be further simplified by performing an Enskog-type 
ordering analysis on all terms to obtain the reduced equation set: [See Appendix, Equation Set 9] 



5]. 



The classical Maxwell thermal velocity twenty moment equation set is: [See Appendix, 



20 



flows. 
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One can extend the above analysis to time average turbulent flows by considering the mass 
mean fluid flow velocity (u) to be composed of a time average velocity (u), and a turbulent motion 
velocity (u'), like Reynolds (1894). The Cartesian components of these velocities are related by: 



The turbulent motion (u') is in viscid, unaffected by collisions since the mean flow velocity 
(u) is conserved in collisions and only the thermal (c) velocities are changed in collisions. This 
important fact has not been previously recognized, explicidy stated or utilized in the 

1 0 scientific/engineering literature. Since both the density and turbulent motion are preserved in 

collisions, there is no coUisional (physical) coupling between turbulent motion fluctuations, density 
fluctuations and thermal velocity correlation fluctuations. Thus, there is no time average correlation 
between these fluctuations. The averaging time required for the time average correlations between 
turbulent velocity fluctuations and molecular (number density and thermal motion correlation) 

1 5 fluctuations to be essentially zero is quite short, since the molecular fluctuations occur on a 

coUisional time scale of nanoseconds and turbulent flow fluctuations occur on much larger fluid 
dynamic time scales (typically measured in milliseconds or greater). Thus, the time dependent time 
average anisotropic turbulent flow equations presented can be used to compute turbulent flows 
changing on fluid dynamic time scales for the vast majority of engineering flows. The time change 

20 adequacy of the time average anisotropic turbulent flow equations needs to be verified for extremely 
rapid flow changes such as those associated with shocks, flames and deteriorations. The only 
turbulent flow variables necessary to analyze anisotropic turbulent flows in the vast majority of 
physical cases are those given in Table I below, all other time average correlations being zero. 



u, = u, +uf 



5 



where now 



V, = Uj +uf +c 
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Table I 
Turbulent Flow Variables 



Huid Dynamic 


Thermal Moments 


Turbulent Moments 


P 


(1) 




(3) 


< 


(3) 




(3) 




(3) 




(3) 






CjCf 


(3) 




(3) 








(6) 




(6) 








(1) 


u;u^u; 


(1) 


Four Fluid Dynamic 
Variables 


Sixteen Thermal Moment 
Variables 


Sixteen Turbulent Moment 
Variables 



The time average turbulent flow variables are the flow density (pj, where the bar indicates 
5 time averaging, the three mass mean velocity components (u^), the three anisotropic thermal kinetic 
energies (c^ , where the two bars indicate both time averaging and averaging over all thermal 
velocities), the three thermal shear components (c,C2 j, the nine thermal kinetic energy fluxes (Cjcf 
and C1C2 ), the thermal structure correlation (C1C2C3 ), the three anisotropic turbulent kinetic energies 
(u[^ ), the three turbulent shear components (u^Uj ), the nine turbulent kinetic energy fluxes (u[u[^ 
1 0 and u[u2^ ) and the turbulent structure correlation ( u[u2U3 ). 

The instantaneous Maxwell molecular velocity (v) and the thermal velocity (c) moment 
equations (Equation Sets 2 and 3) are both valid at every instant of time and the instantaneous 
Maxwell thermal velocity moment equations can be derived from the instantaneous Maxwell 
molecular velocity moment equations. The time average Maxwell molecular velocity (v) and time 
1 5 average thermal velocity (c) moment equations are also both valid simultaneously. The time 

average thermal velocity equations yield equations governing the behavior of the anisotropic time 
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average thermal directional kinetic energies, shears, directional kinetic energy fluxes and structure 
correlations: [See Appendix, Equation Set 10]. 

The time average molecular velocity equations yield equations goveming mass and 
momentum conservation and the behavior of both the time average thermal and turbulent 
5 directional energies, shears, energy fluxes and structure correlations: [See Appendix, Equations 
Set 11]. 

After eliminating the collision change terms with the time average thermal Maxwell 
moment equations and other simplifications, one obtains the equations goveming the behavior of 
the anisotropic time average turbulent directional kinetic energies, shears, directional energy fluxes 
□ 1 0 and structure correlations: [See Appendix, Equation Set 12]. 

Ln Both sets of third order directional kinetic energy fluxes and structure correlation equations 

l^y are closed using the fourth order density gradient independent time average thermal and turbulent 

V moment closure relationships: [See Appendix, Equation Sets 13 and 14]. 

''S One thus obtains the complete closed Maxwell twenty moment equation set goveming 

I^f 1 5 anisotropic time average turbulent flows: [See Appendix, Equation Set 15]. 

P By the same process, corresponding higher order time average turbulent moment equation 

sets may be derived using higher order closure relationships. The general density gradient 
independent time average thermal and turbulent even moment closure relationships are: [See 
Appendix, Equation Sets 16 and 17]. 

20 These closures allow the infinite set of time average Maxwell moment equations goveming 

turbulent flows to be closed at any desired odd order level starting at the third order. If one closes 
at the third order (say), then one can solve all higher order moment sets (fourth and fifth, six and 
seventh, etc.) in order, where the errors are still of Burnett order in the convective terms of the 
highest order moment equations (fifth, seventh, etc.) considered in the set. The method thus allows 
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solution of all the Maxwell moment equations, both those sets below and above the chosen closure 
level. 

The above results may be simply extended in many ways. For example, the general even 
order density gradient independent moment closure relationships (Equation Set 4) are valid in 
5 equilibrium. Thus, [See Appendix, Equation Set 18] and the even order coUision integrals 



It will be recognized that the above-described invention may be embodied in other specific 
forms without departing from the spirit or essential characteristics of the disclosure. Specifically, 
the anisotropic turbulent flow equations may be reformatted and/or approximated and computed in 

1 0 many different ways to solve specific problems. Obvious first steps would be to add gravitational 
forces or to sum the directional turbulent energy equations to obtain a total turbulent energy 
equation or to reduce the complexity of the equation set through a Prandtl (1904) boundary layer 
analysis. The analysis presented (valid for monatomic perfect gases) can be expanded to diatomic 
(and polyatomic) gases, gas mixtures, reacting gas mixtures, plasmas, etc. through known methods, 

1 5 Chapman and Cowling (1970). The analysis can be extended to isotropic fluids by setting the three 
directional thermal kinetic energies cf or c^ equal, and solving only the summed total thermal 
energy equation, rather than the three directional thermal kinetic equations. Since this reduces the 
Equation Set 9 to the Navier Stokes equations, Kliegel (1990), the analysis can be extended to 
liquids by setting both the three directional thermal kinetic energies equal and the density constant. 

20 The closed anisotropic turbulent flow equation set can be utilized with existing Reynolds averaged 
Navier Stokes equation solvers to solve incompressible isotropic fluid problems. Thus, it is 
understood that the invention is not to be limited by the foregoing illustrative details, but rather is to 
be defined by the appended claims. 




Cj C2 C3 - Ci C2 C3,etc. may be calculated for all molecular models. 
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While all equation sets may be solved by straightforward computation, they may be more 
efficiently solved using the method described in my co-pending application entitled "Method for 
Algebraically Solving Differential Equations, Including Stiff Equations, to High Accuracy". 



FILE NO. 004210.P001 



-17- 



GWHVDLM 



APPENDIX 



Equation Set 1 
(1-1) 



^[pF(v)] = l[pF(v)].v.A[pp(,)],.^^[pP(.)],V3A[pp(,)] 



pi;[F(v)*-F(v)] 
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Equation Set 2 
(2-1) 



|[p v;v-v^]-h^[p v;xv^]+^[p v[vr<]^^[p v!v-vr^ 



9k. 



= P^[v!vX*-v|vxj 



where vjvj V3 is the instantaneous molecular velocity moment vjvj V3 averaged over F(v) and 

* * 

vjv2 V3 is the instantaneous scattered molecular velocity moment vjv2 V3 averaged over F(v) . 
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Equation Set 3 
(3-1) 



|-[c?] + £hc,] + ^[c.C3] 



ok 



£[c,C3]-H£hc3]-f^[c3^] 



+ / 



da 



da, 



^1 *^2*^3 -V ^*"1^2 *"3 -V ^*^1^2*^3 n 



^r::z; 



c, 0203- 



/(cf) C,''-'C-C3"-^[C,C,] c;crC3'-/7[c.C3] c^C-C,"-^]!^ 

Cfc-C3"-^1-^ 



+ [cfcrc3"^'-/[c,C3] c';'c1cl-m\c^c,\ c[cr'<^l-nY\ cfc^c 



n-J 
3 



p 

1 ap 



p 8x3 



where 0(0203 is the instantaneous thermal velocity moment cjcj C3 averaged over F(v) , and 

* if. 

c,C2 C3 is the instantaneous scattered thermal velocity moment CjC^Cg averaged over F(v) . The 

latter may be calculated from the full Boltzmann collision change integrals, Kliegel (1990). 



FILE NO. 004210.P001 



-20- 



GWH\DLM 



Equation Set 4 
(4-1); (4-2); (4-3) 



m n-1 
'2^3 



neglected terms being of Burnett order, where /H-m+n = 2N +l,an odd number. 
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Equation Set 5 
(5-1) 



+ / 



+ m 



CjCjCj 



+ c; C2 C3 



-^[crV^c^] + c,c, -^[c|-^c^"cj] + c.C3 ;|-[c;-V;c;] 



<^3 dx.^ 



^1 ^2 ^3 :x 



+ CJC2C3 +CjC2 C3 



+ n c. C7C3 



* 1 

C1C2C3 -C,C2C3 



^+c;cr'c3 



^ ^3 



«3x3 



+ C,C2 



<3x 



^0 



where / + m+n = 2N +l,an odd number and the neglected terms are of Burnett order in the 



convective terms (left hand side). 
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Equation Set 6 
(Twenty Moment Equations) 



Continuity (6-1) 



Momentum (6-2) 



Directional Thermal Energy (6-3) 




Thermal Shear (6-4) 
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Equation Set 6 

Directional Thermal Energy Fluxes (6-5) 



# 
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(6-6) 



Equation Set 6 



3 r — 2 



+ u, 



Bx, 



C.C^] + U3-^[c,C^] + U3-i-[c,C^] 



+ C,C2 



— 5" 3u 



3x, 



+ C2C2 



dx. 



■ + C3C2 



9x, 



+ 2 I c,c? — ^ -H c,c, — ^ -I- c,c,c 



axjj 



- 2 CjCj 



8x 



dx. 



:j r::2i o — r— ii i 

pax, 



[cfc^-cf [c^]-2c,C2 [c.Cj]] 



+ |C,C2C2 -3 C1C2 



+ 



p dxj 



m 



[c,C3C^-C,C3 [c^]-2c,C2 [C2C3]] 



— 11 1 dp 

p dxj 
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Equation Set 6 



Thermal Structure (6-7) 



9x, 



+ C1C2C3 



duj 



H-CX 



3^2 



+ C.C 



3^1 



2 

ax. 



+ C,C2C3 



8X2 

9ui ^ 8U| 
8x2 ^ ^ 3x3 
9u, — T du 



8x2 ^ ^ 8x3 



2^1 



- C2C3 



2^ + ^^ 

3xi ^ ^ 9x2 
3 r^i a 



+ C1C2C3 



3u3 
8X3 

a r 1 



ax 



C1C3 



- C,C2 



a r^i a 



3-hc2]-^^[c^]-h— [C2C3] 



dx 



dx. 



ax. 



+ _|_[c^C3cf] + ^[c,C3C^] + A[c,C3C^] 



ax 



3x, 



[c,c3c^ -e;^; [c^]-2 s;^; [c^]]^^ 
[cjCjc' - ^ [c^] - 2 e;^; [3^]]-^!^ 



= C,C2C3 



where =cf + Cj + Cj and the collision change terms (right hand sides) have been calculated 
using the method of Kliegel (1990). 
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Equation Set 7 
(7-1; 7-2; 7-3; 7-4) 



c.Cjcf = 3 cjc; [cf ] 

cfc2=cf [c^]+2c;c;[c;c;] 
C1C2C3 = c;c^ k 1 -I- 2 cjc^ [c^] 
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Equation Set 8 
Closed Directional Themal Energy Fluxes (8-1; 8-2) 
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Equation Set 8 



Closed Thermal Structure (8-3) 




3 

the neglected terms are of Burnett order in the convective terms (left hand sides). 
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Equation Set 9 



Continuity (9-1) 



^+3^[puJ-h3^[pu,]-h^[pu3] = 0 



ax 



Momentum (9-2) 



du, . 9u 

"at 



9u, 



au, 



- + Uj — + U2 ~- + U^ - 



axj ^ ax2 ^ axg 



1 

+ — 
p 

= 0 



Directional Thermal Energy (9-3) 

a r^i . a ["Ti . a r"? 



1 :s +^i'^2 ^ 
dx, dXj 



dx. 



t _ ^ _ ^ _. 



+ — 
P 



Thermal Shear (9-4) 



c,C2 = — 



V 



■+c: 



dx. 



Directional Thermal Energy Fluxes (9-5) 



c,c, = - 



2v 



I- J ax, 



2 ^2 

6 1; 



Equation Set 9 



axo 
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Thermal Structure (9-6) 

neglected terms being of Burnett order in the thermal shear, directional energy fluxes and structure. 
The shear stress and thermal energy flux relationships are identical to the results Kliegel (1990) 
obtained by an Enskog-type solution of the Boltzmann equation. 
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Equation Set 10 
Time Average Directional Thermal Energy (10-1) 

at 









- a 




Hl^^'ax. 


h'l^^^ax. 




+U3— 


of' 








0X3 



+ 2 



=j au, 



au 



au, 



dXi 0X3 



P 



cf 1 + 



ax. 



ax 



a r-==f 

— [p C3C, 



Time Average Thermal Shear (10-2) 



+ c,c 



— =f 9u, 



9u, 

* ^ 3x, ^ 8x2 ^ ^ 8x3 



= 8u2 



+ 3-^ + CiC2 — 
dx, dXj 



3x. 



1 


" a 




^ a 






+ = 

p 


ax, 




axj 




+ A[pe,c,C3] 



1^2 
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Equation Set 10 



Time Average Directional Thermal Energy Fluxes (10-3) 
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Equation Set 10 

Time Average Directional Thermal Energy Fluxes (10-4) 
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Equation Set 10 



Time Average Thermal Structure (10-5) 
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Equation Set 1 1 



Continuity (11-1) 



Momentum (11-2) 



+ ^[p[uiU2+u;u^+c,C2]] 



= 0 



Turbulent Directional Energy (11-3) 



— n In. +v'^ . f ' 

at I 



-iPiur-Hur+cf]]+A 



+ p jugU/ +U3U[^+2UiU[U3+U3U[^ +U3C^ +2UiC,C3 +036^ 



= pv 
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Equation Set 1 1 



Turbulent Shear (11-4) 



^j^p|^U,U2+u[U2 +C1C2 ]]'*"^~[^ [^2^1^ ■•"^2^1^ "*"^2^[^ H-2u,u[U2 + Uj + 2 Uj C,C2 jj 

^ T— r 2 — "72 , o,, „/,/ . , /,/2 . . rrU . 2 



3x 



P U1U2 +U,U2 +2U2U|U2 + U[U2 + 2 U2C,C2 + U1C2 +C,C2 



2 

+ [uiU2U3+U3U[U2 + U2U[U3-|-UjU2U3 + u[U2U3 +U3CjC2+ U2CJC3 -hUjCjCj + C1C2C3 jj 



1] 



9x^ 

= - pD c,C2 

Turbulent Directional Energy Fluxes (11-5) 



^1 P 1^1^ UiU[^ +u|u[^+ 3Uj cf 4-CiC 



+^[p [^^;' +6 {^;'u;' +4 n^u[u{' + u{\' ]] 

]] 



^ r~ Itt"^ . ^ rr2„/2 . A 7r„/„/2 . „/2 

^ T-f^— 2=2 . ^ ^2 . ^„/2^2 . ^2^2 



4-— I p|6ui c, +4uiCiCj +6ui Cj +Cic; 

OX I 



^^l^P [^Ui^U2 +3Ui^u[U2 +3UiU2U[^+U2U[u[^ +3Ui U2U[^ + U[U2U[^ jj 



a r-r.- 



+- — p 3ui c,C2 +3 UjU2Cj -i-U2C,Ci + 3U1C2C1 H-3u[ C,C2 +3u[u2Ci +C1C2C 
dX2 L I- 

+ -^[p |u^^u^4-3u^^u[u3 +3u^u^u[^ -hu^u(u[^ +3u^uX^ + u[u3u[^ jj 



+ - |p|3Ui CjC3 + 3 U1U3C1 +U3CiCj +3U1C3C1 +3u[ C(C3 +3u[U3Cj +C1C3C 



1] 



= pv 



3u. 



1 

+- 
2 



c,c -3 c,c, 
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Turbulent Directional EnergTRuxes (11-6) 



p |^UiU2"+U,U2^ -|-2U2U[U2+U[U2^ '^^\^2 + 2 U2C1C2 -hCjCj jj + 

3 r-r— 2— 2. —2-;?. 2—^ 



9_r_r 2 

at I 



-h [ Ui^U2^ + U,^U2^ -|-U2^U[^ -|-4UiU2U[U2 "h u[U2^ + 2U2 ^U^^ + uf jj 



3 f— r — 2T — 2 "T 



+ -^1 p| Ui"C2+U2"cf + 4UiU2C,C2+2u,CjC2+u[^C2+2U2C2Cj' + U2''c;' + 4u[U2CiC2 + <^{^2 



,'2^2 



.2„2 



1] 



+^[p [u.Uj' +3u2^u;u2 -i-3u,U2U2^-i-u,uX^ +3u2u;u2^ + u;uX^ ]] 



d r-r.— 2 



]] 



-h-^lpl 3U2"C,C2 +3 U1U2C2-I-U1C2C2+ 3U2C1C2 +3U2^CiC2 -|-3u[U2C2 +C,C2C2 

8x2 L L 

H [ UiU3U2^+U|U3U2^-|- Uj^U^Uj +2 U1U2U2U3 +2 U2U3U[U2 -|-2U2U[U2U3 + UgUj^ + U3 u[U2^ + u[U3U2^ 

3x. L L 



H-^ \p\vi2 CiC3-|-U,U3C2-|-2UiU2C2C3 +2U2U3C1C2 +UjC3C2 + 2U2C1C2C3 + 



U3C,C^]j 



a r-i 



H-^ 1 p| U2 CjC3 +2u(U2C2C3 + u[U3C2 +2U2U3C,C2 + C1C3C2 



1] 



= pv 

















3 ' 



1 



-lUjCjCj +^|c,C^-9 C|C2 
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Turbulent Structure (11-7) 

[^UjUjUg H-U3C,C2 + U2CiC3 +U1C2C3 4* C,C2C3 + U3 U^Uj + U2U[U3 + U,U2U3 + u(U2U3 Jj 
+ ^^|p I U2U3u/-|-Ui^U2U3 +2 UjU2U[U3 +2 U,U3U[U2 +U2U3U[^ +2UjU[U2U3 +U2 UjU^^ + U3 UjU^^ jj 

+ ^^j^p|u2U3U[^ 4- Ui^C2C3 +2u,U3C,C2 +U2U3cf +2U1U2C1C3 +2UiCjC2C3 +U2C3cf + U3 C2cf jj 

+ ^^j^p|u[^C2C3 +2u[U2CjC3H- 2u[U3CiC2 + UjUjCf +02030^ jj 

H-^^|p 1^ UiU3U2^+U2^u(U3 +2 UjUjUj^ +2 U2U3U[U2 +U,U3U2^ +2U2U[U2U3 +U, U3U2^ + U3 u[U2^ jj 
H-^^j^pj^u[U3U2^ + U2^CjC3 -I-2U2U3C1C2 +U,U3C2 + 2Ui U2C2C3 +2U2C1C2C3 H-Uj C3C2 + U3 C1C2 jj 
+ ^^|^p|uf C1C3 +2u[U2C2C3+ 2U2U3C1C2 + U[U3C2 +CiC3C2jj 



+ -^|p ^ UjU2U3^+U3^u[U2 +2 U2U3U[U3 +2 U^UjUjUj -|-UjU2U3^ -|-2U3U[U2U3 +U2U[U3^ + U, U2U3^jj 



dx 

d 



+ ^ |^p|ui^2^3^ + ^3^^1^2 +2U,U3C2C3 H-UjUjCj +2U2U3C1C3 4-2U3C,C2C3 H-U2CiC3 + U, C2C3 jj 



j^p^uf C,C2 +2U2U3C,C3+ 2u[U3C2C3 + u[U2C3 +CiC2C3jj 

= -pvj^u 



1 + ^2 C1C3 + U3 C1C2 + -C1C2C3 
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Equation Set 12 



Directional Turbulent Energy (12-1) 




Turbulent Shear (12-2) 
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Equation Set 12 
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Equation Set 12 
Directional Turbulent Energy Fluxes (12-4) 



- a 



u;uf 1+U3 



—JL 

8x, 



u;uf] 



+ 2 



' ' ' ' dx, ' ' ' 8X3 



ox, dXj 



- 2u;u2 



3x, 



- u. 



8x 



3x, 



dx 



8x, 



+ 2 



ax 

a 



7777 a 



ax 



ax. 



ax. 



7~7 a [ 2! , „/2 a r 2I , , /,./ a \ 2 



+ U,U2 l + U, 



dx. 



^x^ 



C2 +u,U3 - — icj" 



ax. 



= 0 
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^[u;u^u;] + u, ^[u;u^u;] + ^[u;u;u;] + A[u;u^u^] 



-7-7-7 



du 



+ U;U2U3 — + UjUj^ — 

dx, dx 



2 

du 



/,,/2 



+ U3U1 



^ . „/„/„/ 



+ U1U2U3 



dx 



8x3 

OX. 



, ,2 8U 



dx. 



— U2U3 



dx, 



dx 



7-7 d 



9xj 



3x 



3x. 



3xj 



l + CtC3 



+ C,C3 



3x, 



+C2C3 



3x2 ^ ■' 3X3 



[u;u;] 



+ c 



~2 a r / / I ^= a 



ax 



ax. 



= 0 
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Equation Set 12 TuAulent Stnicture (12-5) 



GWH\DLM 



Equation Set 13 (13-1; 13-2n3-3; 13-4) 



c?c2=cf [c^]+2c;c;[c;c;] 

c.CjC^ = c^2 [cs] + 2 c;c; [c^c^] 
neglected terms being of Burnett order. 

Equation Set 14 (14-1; 14-2; 14-3; 14-4) 



neglected terms being of Burnett order. 
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Equation Set 15 
Time Average Continuity (15-1) 

Time Average Momentum (15-2) 

a r- 



+ ^[p[u^u;+u;u;-i-s;^]] 



"■3 
= 0 
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Equation Set 15 
Time Average Directional Thermal Energy (15-3) 
d_ 

at 







— a 























+ 2 



+ CtC 



1^3 



9x. 



Time Average Thermal Shear (15-4) 



OX, ox, ox. 



+ C1C2 



=T au 



+ cr 



1 

+ = 
p 



au. 



au. 



ax 



dX2 0X3 



ax. 



ax. 



pc.Cj 



+ 



a r- = 

a;r^^' 



C2C3 
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Equation Set 15 

Closed Time Average Directional Thermal Energy Fluxes (15-5) 




FmENO. 004210.P001 
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Closed Time Average Thermal Structure (15-7) 




FILE NO. 004210.P001 



Equation Set 15 



Time Average Directional Turbulent Energy (15-8) 



+ 2 



— du, —y-jdu. —TT^^i 
U, rr-^ + uJUj r— ^ + u'u{ - 



M"3 



3x. 



1 
P 

= 0 



^ r-„/„/2i 



9xj 



3x 



Time Average Turbulent Shear (15-9) 



-h U1U2 



7-7 au. 



3xi 



+ u 



/2 



au 

3x 



1 —r-j^^] 



-h u 



/2 



9u2 
3^ 



-h U1U2 



axn 



-i-u;u3 



3x, 



1 
p 



9x2 - 



Bx. 



pu[u;u;] 



= 0 
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Equation Set 15 

Closed Time Average Turbulent Directional Energy Fluxes (15-10) 




FmENO. 00421 0.POOl 
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Equation Set 15 

Closed Time Average Turbulent Directional Energy Fluxes (15-11) 
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Equation Set 15 



t n 



Closed Time Average Turbulent Structure (15-12) 

+ u u^u ^ + T-^ + u^u^' — L 
dx, dXj 0X3 

' ' 8x, ' ' ' 9x2 ' ' 9x3 

. . ^7-7-7 



.'..'2 



dx, dXj dX3 

= 0 
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Equation Set 16 



(16-1; 16-2; 16-3) 



neglected terms being of Burnett order. 



Equation Set 17 

(17-1; 17-2; 17-3) 

ufXV='P] u;'-'u>r+'«[<u^] w'u^^^v+^lulu;] uf^v 

u;V*V=4^1 u;'-'u^'"u;" + m [uf ] u;'u^'"-y"-i-n[u^] ufuf< 

U, U2 U3 =/[UiU3j Ui U2 U3 +m[U2U3j Ui U2 U3 + H j^Uj j Uj U2 U3 

neglected terms being of Burnett order. 
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Equation 18 
(18-1; 18-2; 18-3) 



# 



^1 *^2^3 ~ 



C,C2 C3 



C1C2C3 



/ 

= 0 
= m 
= 0 
= n 
= 0 



c{ ^CjCj, / odd, m and n even 



/, m and n all odd 



C1C2 C3, m odd, / and n even 
/, m and n all odd 



c;c2C3 , n odd, / and m even 
/, m and n all odd 



neglected terms being of Burnett order, where / + m+n = 2N +l,an odd number. 
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